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We investigate the dynamical properties of a class of spatially homoge- 
' neous and isotropic cosmological models containing a barotropic perfect fluid 

and multiple scalar fields with independent exponential potentials. We show 
^P' that the assisted inflationary scaling solution is the global late-time attrac- 

ts ■ tor for the parameter values for which the model is inflationary, even when 

d ' curvature and barotropic matter are included. For all other parameter values 

the multi-field curvature scaling solution is the global late-time attractor (in 
these asymptotically stable solutions the curvature is not dynamically negligi- 
ble) . Consequently, we find that in general all of the scalar fields in multi- field 
models with exponential potentials are non-negligible in late-time behaviour, 
contrary to what is commonly believed. The early-time and intermediate be- 
haviour of the models is also studied. In particular, n-scalar field models are 
investigated and the structure of the saddle equilibrium points corresponding 
to inflationary m— field scaling solutions and non-inflationary m— field mat- 
ter scaling solutions are also studied (where m < n), leading to interesting 
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transient dynamical behaviour with new physical scenarios of potential im- 
portance. 

PACS numbers(s): 04.20. Jb, 98.80.Hw 
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I. INTRODUCTION 



Inflation is generally considered to be a reasonable solution to many of the fundamental 
problems within the standard cosmological model. In original proposals the early uni- 
verse experiences a period of accelerated expansion and essentially expands at an exponential 
rate, (i.e., R{t) oc e"^* where R(t) can be considered as the size of the universe and A is some 
positive constant). Since these early stages, there have been a variety of inflationary models 
that include scalar fields that have been proposed 0, and scalar fields have come to play 
an important role in determining the dynamics of the early universe. In one important class 
of inflationary models the condition of exponential expansion is relaxed, and the universe 
grows at a power-law rate, R{t) oc t^, where p > 1 In particular, power-law inflationary 
models arise in models with a scalar field having an exponential potential V{(j)) = Voe'"^ 
1^. Although power-law inflation is successful in solving the horizon and flatness problems, 
inflation in these models persists into the indefinite future and a phase transition is required 
to bring inflation to an end (however, see [^]). 

Spatially homogeneous models containing a scalar field (p with an exponential potential 
have been analyzed extensively 0,0]. It is known that all ever-expanding scalar field mod- 
els experience power-law inflation when the parameter /c^ < 2; i.e., when the potential is 
sufficiently flat. The models have also been studied when k"^ > 2 [0. Recently cosmological 
models containing both a scalar field with an exponential potential and a perfect fluid with a 
linear barotropic equation of state have been studied. It is found that in the general class of 
Bianchi type B models that the power-law inflationary solution is still the global attractor in 
the physically realistic regime (i.e., when 7 > 2/3) if fc^ < 2 Interestingly, the addition 
of a barotropic perfect fluid creates the existence of a new type of solution, appropriately 
called a matter scaling solution PHTTD , in which the energy-density of the scalar field scales 
with that of the matter; the effective equation of state for the scalar field is the same as that 
of the perfect fluid. The stability of the matter scaling solution has been studied in P,p^. 

Exponential potentials arise in many theories of the fundamental interactions including 
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superstring and higher- dimensional theories P,[T2[]. Typically, 'realistic' supergravity theories 



predict steep exponential potentials |jT2|] (i.e., /c^ > 2), effectively eliminating the possibility 
of power-law inflation. However, dimensionally reduced higher-dimensional theories also 
predict numerous scalar fields, and so it is of interest to study models with multiple scalar 
fields. 

In the recent work of Liddle, Mazumdar and Schunck [l^ the effect of additional scalar 



fields with independent exponential potentials was considered. They assumed n scalar fields 
in a spatially fiat Friedmann-Robertson- Walker (FRW) universe. They found that an arbi- 
trary number of scalar fields with exponential potentials evolve towards a novel inflationary 
scaling solution, which they termed assisted inflation, in which all of the scalar fields scale 
with one another (and are hence non-negligible asymptotically) with the result that infla- 
tion occurs even if each of the individual potentials is too steep to support inflation on its 
own. The existence of multiple uncoupled scalar fields, each having an exponential potential, 
could therefore, through a combined (or assisted) effort, be a source for power-law inflation. 
This is true even though each individual scalar field need not be a source for inflation, and 
might therefore lead to compatibility with supergravity theory. 

In a recent dynamical analysis it was shown that this assisted inflationary solution is 
a late-time attractor in the class of zero-curvature FRW models. This was done by choosing 
a redefinition of the fields (a rotation in field space) which allows the effective potential 
for field variations orthogonal to this solution to be written down; in analogy with models 



of hybrid inflation [|T^ it was then shown that this potential has a global minimum along 
the attractor solution. Also, analytic solutions describing homogeneous and inhomogeneous 
perturbations about the attractor solution without resorting to slow-roll approximations 



were presented in and curvature and isocurvature perturbation spectra produced from 
vacuum fluctuations during assisted inflation were discussed. 

In this paper we shall present a qualitative analysis of models with the action 



n n 
i=l i=l 



+ Sm, (1.1) 



where Sm is the matter contribution. Almost all previous analyses of multiple scalar field 
infiationary models have assumed zero- curvature FRW spacetimes with no matter; here we 
extend the analysis to include both curvature and matter. In section ^ we shall present the 
governing equations for n scalar fields with exponential potentials and matter. In section 
pn] we shall study the two-scalar field model with no matter, and, in particular, discuss 
the stability of the two-field assisted infiationary model. In section |^ we shall study the 
two-scalar field model with barotropic matter. In section we shall discuss three- and 
multi-scalar field models. In section ^ we present our conclusions. 



A. Generalized Assisted Inflation 

Recently, models with n x m scalar fields (pij and containing multiplicative exponential 
terms in the effective potential of the form 

n m n 
1=1 j = l i=l 

where 1 < i < n and 1 < j < m and kij are n x m real positive constants which are not 
zero, have also been studied. A qualitative analysis of the m = 1 case has been given in 
, where an analogy was made with the dynamics of soft inflation . 



In a class of spatially flat FRW multi-scalar flelds models with multiplicative expo- 
nential potentials were studied. Potentials of this form are quite common in dimensionally 



reduced supergravity models [19,20]. Exact two-fleld and general n-fleld power-law scaling 



inflationary solutions were obtained, which were demonstrated to be late-time attractors. 



generalizing the assisted inflationary solutions previously obtained [|T3|; this behaviour was 
dubbed 'generalized assisted inflation'. It was shown that it is more difficult to obtain as- 
sisted inflation in these generalized models with cross-couplings between the scalar flelds in 
the potential; the flelds in any one exponential term tend to conspire to act against one 
another rather than assist each other (a result also noticed in EOl). However, as with the 



original version of assisted inflation, this inhibiting affect can be compensated for if there 



are enough exponential terms present in the potential (i.e. if n is large enough) [1T3| , |T8 



The dynamics of 'generalized assisted inflation' was investigated in more detail in pT 



By introducing field rotations, which results in the introduction of two orthogonal fields one 
of which is massless and the other posseses an exponential potential [0,^], the nature of the 
late-time attractor solution in a particular class of models was determined. A dimensionally 
reduced action resulting from within the context of a generalised toroidal compactification 
of higher-dimensional fields in Einstein gravity minimally coupled to massless scalar fields 
was shown to give rise to a model of the form under investigation, and it was shown how 
the addition of interactions between the fields impede inflation in this model. 

Similar behaviour was also noted by Kanti and Olive [^] in multi-field assisted inflation- 
ary models with standard chaotic polynomial (rather than exponential) potentials, which can 
arise in modern Kaluza-Klein theories (and are a natural outcome of the compactification 
of higher dimensional theories down to four dimensions). Indeed, Kanti and Olive ||2^ have 
recently proposed a possible realisation of assisted inflation based on the compactification 
of a five-dimensional Kaluza Klein model, and have shown how the additional fields of the 
assisted sector actually impede inflation (they also showed that the assisted sector, coming 
from a Kaluza Klein compactification, eliminates the need for a fine-tuned quartic coupling 
to drive chaotic inflation). In Kaloper and Liddle [^| the dynamics of a simple implementa- 
tion of the idea in Kanti and Olive |^ was analyzed in more detail. Since assisted inflation 
no longer corresponds to an asymptotic attractor, they found that as inflation proceeds the 
number of fields participating in the assisted behaviour decreases resulting in the interesting 
novel feature that the density perturbations generated retain some information about the 
initial conditions. 



II. THE MODEL 



We shall assume that the spacetime is spatially homogeneous and isotropic. The line 
element for such a spacetime has the form 



-dt' + R\t) 



+ r\de' + siYi\e)dr 
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where k = +1,-1,0 determines whether the model is closed (positive-curvature), open 
(negative-curvature), or flat (zero-curvature). 

We shall consider n scalar fields (pi, where 1 < i < n, in which the effective potential has 
the form 



1=1 



where the fcj are real non-zero positive constants. We also assume that there exists a non- 
interacting perfect fiuid with density p and pressure 

P={l-l)p, (2.1) 

and we shall assume that 1 < 7 < 2. The Einstein field equations, the conservation equa- 
tions, together with the Klein-Gordon equations for the scalar fields, yield the following 
autonomous system of ordinary differential equations: 



H = -H'-l(^J2 + l^^ff - ^(37 - 2)p, (2.2) 

p = -3-fHp, 

where = k/B? is the curvature of the spacelike hypersurfaces, H = R/R is the Hubble 
expansion, and an overdot represents differentiation with respect to coordinate time t. Units 
have been chosen so that StcG = c = 1. 

To analyze the dynamical system given in equation (|2.2|) we transform the system using 
expansion normalized variables. Expansion normalized variables have had great success in 
analyzing the asymptotic behaviour of many cosmological models. See for arguments in 



support of using dimensionless expansion normalized variables. One primary reason involves 
the decoupling of one of the differential equations, which effectively reduces the dimension of 
the system by one, and, in some cases, compactifies the phase space. We choose expansion 
normalized variables of the form 



The resulting dynamical system describing these perfect fluid multiple scalar field models 
becomes 

— = ^](2g-37 + 2), 

'^^%{,-2)-^K^l (2.4) 
for (1 < « < n), where the deceleration parameter has the following form 



i=l i=l 

and 

n n 



i=l i=l 

Assuming a non-negative energy density (i.e., f2 > 0) and if < 0, (i.e., in the negative 
and zero-curvature cases) the phase space for the dynamical system in the expansion nor- 
malized variables (f2, $j, ^>i) is compact. If ^i? > (i.e., in the positive curvature case) then 
the transformation given by equation (|2.3|) becomes singular when H = 0. Here we shall 



only make some partial comments with regards to the asymptotic behaviour of the posi- 
tive curvature models. All of the equilibrium points correspond to self-similar cosmological 
models and hence to power- law solutions 

III. QUALITATIVE ANALYSIS OF TWO-SCALAR FIELD MODEL 

We shall first discuss the dynamics of the model with only two minimally coupled scalar 
fields and with no matter. We obtain this model by setting n = 2 and f2 = in ( |2.4| ). In 
this case we obtain the four- dimensional dynamical system given by: 



8 



where 



and 



^^2 T / ON V^, ^ 2 
-— = ^2(g-2) - — /C2$2 
CLT Z 

dr 2 

^ = ^2{q+l + ^h^2) (3.1) 



g = 2^1 + 2^2 - - $2 



3^ 

-1 + ^1 + ^2 + "^l + ^l- 



It is possible to choose simphfied variables as in pH] via a rotation in field space; although 



this would simplify the analysis of the assisted infiationary solution, it would perhaps be 
more difficult to discribe all of the qualitative properties of the models and relate this analysis 
to previous work. 



A. Assisted Inflation 



The fiat Assisted Inflation model corresponds to the equilibrium point A of the 
system (|3.1| ) given by 

{^l,*2,$l,$2}^ = 

; kih' h^h , v/6(fci^ + k,^) - k,^k,^ v/6(fci^ + fc2^)-fciV ^ 



' ^/Qik^' + k^^y V6{k^^ + k^^y V6{k^^ + k^^) ' V6{k^^ + k^ 

(3.2) 

which is equivalent to 



\/6ki V6k2 \/6ki VQk2 
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where 

The deceleration parameter for this solution is given by 



(3.3) 



and hence this solution, with 



R{t) oc 



and 



is inflationary (g^ < 0) if 



2 

P 



2 V ^ = 2K-^ = — i — > 1; 2 > K^. (3.4) 

Since a single scalar fleld can only give rise to an inflationary power-law solution if 
^ > i for i = 1 or 2 |,|], this means that the two-scalar fleld model can be inflationary 
even when the each of the individual potentials is too steep for the corresponding single 
scalar fleld model to inflate (and hence the terminology assisted inflation). The eigenvalues 
corresponding to the equilibrium point A are given by (See Appendix for details.) 

- ({K^ - 6) ± y^iK^ - 6)2 + SK^iK^ - 6)) (3.5) 
2 4 \ / 



Hence this equilibrium point is stable when (|3.4| ) is satisfled, and so the corresponding 
assisted inflationary solution is a late-time attractor . 
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B. Stability of Equilibria 



We note that several of the equihbrium points occur in the three-dimensional invariant 
set corresponding to the zero- curvature models defined by 

When matter is included, there exists a monotonic function so that in the full dynamical 
phase space there can be no periodic or recurrent orbits and the global dynamics can be 
determined. This implies that the qualitative features described in this section can be more 
rigorously proven. All of the equilibrium points and their corresponding eigenvalues are 
listed in Table |. Using this this table let us discuss the local stability of these equilibrium 
points. 

As noted above the equilibrium point A, given by (|3.2|) , corresponds to the assisted 
inflationary solution. It exists for all parameter values satisfying 

^ + '^'^^ 

and is a sink (late-time attractor) for all parameter values satisfying ( p.4|) (else it is a saddle). 

There are two equilibrium points, denoted by Pi and P2, whose coordinate values and 
associated eigenvalues are given in Table |, which correspond to zero- curvature power-law 
solutions in which one scalar fleld (either 0i or 02, respectively) is negligible; these solutions 
exist if I < and are inflationary if, in addition, | < (for each i = 1,2, respectively) 

i i 

and correspond to the well-known single scalar fleld power-law solutions . From Table | 
we see that each Pj has two negative eigenvalues and one positive eigenvalue for all relevent 
parameter values and an additional eigenvalue which is negative if kf < 2 (and positive for 
2 < kf < 6); hence these points are saddles and have a one- or two-dimensional unstable 
manifold depending upon whether kf < 2 ot kf > 2, respectively. 

There also exist equilibrium points, denoted by C^i, CS2 and CS, whose coordinate 
values and the associated eigenvalues are given in Table |I[ The solutions correspond to 
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power-law solutions in which the curvature scales with the first scalar field, the second 
scalar field or both, respectively. The single- field curvature scaling equilibrium points CSi 
and CS2 are both saddles. The two-field curvature scaling equilibrium point CS is a sink 
whenever | > + ^ (otherwise a saddle) . Whenever the two-field curvature scaling solution 
is stable, it necessarily has negative curvature. 

There is an equilibrium point, denoted by M, corresponding to the Milne form of fiat 
spacetime, which is always a saddle. 

Finally, there is a one-dimensional set of equilibrium points parametrized by \l/o, denoted 
by MSF, corresponding to zero-curvature massless scalar field models (in which both po- 
tentials are zero). There is one zero eigenvalue corresponding to the fact that there is a 
one-dimensional set of equilibrium points. There are values for \l/o for which the remaining 
three eigenvalues are positive and hence a subset of MSF are sources (the remainder are 
saddles). These correspond to well-known early-time attracting massless scalar field models 
i- 

C. Discussion 

From the analysis above we conclude that the two-field assisted infiationary solution A 
is the global attractor when XlLi > k ^"^^ two-field curvature scaling solution CS 
is the global attractor when X]^=i < |- "^^^ massless scalar field solutions MSF are 
always the early-time attractors. 

In all cases both scalar fields are non- negligible in generic late-time behaviour. This is 
contrary to the commonly held belief that in multi-field models with exponential potentials 
the scalar field with the shallowest potential (i.e., smallest value of k) would dominate at 
late times. Indeed, we have shown that the single field power-law infiationary models always 
correspond to saddles, so that we have the rather surprising result that generically a single 
scalar field model never dominates at late-times! 

We note that both the assisted inflationary solution and the massless scalar field early- 
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time attractors correspond to zero-curvature models. However, the curvature is not always 
dynamically negligible asymptotically because the two-field curvature scaling solution has 
non-zero curvature. 

There is a range of parameter values for which the assisted infiationary solution is the 
global late-time attractor (when the solution is non-inflationary it corresponds to a saddle). 
For all of these parameter values the single field power-law solutions Pi and P2 are sad- 
dles. However, there are allowable parameter values for which either Pi and P2 are both 
infiationary, or one is infiationary while the other is not, or and both are non- infiationary. 
This might give rise to some new interesting physical scenarios. For example, a model could 
asymptote towards an infiationary single field solution Pj, stay close to Pj for an arbitrar- 
ily long period of time (since Pj is an equilibrium point) infiating all the time, and then 
eventually leave Pj and evolve towards the stable attracting infiationary solution A. (Note 
that if either of Pi or P2 are infiationary, then A is necessarily infiationary -see p.4|) . This 
is akin to a double-infiationary model |]2^ in which the density fiuctuations on large and 



small scales decouple (i.e., the scale invariance of the spectrum is broken) thereby allowing 
the possibility of more on large scales which is in better accord with observations. 

IV. QUALITATIVE ANALYSIS OF TWO-SCALAR FIELD MODEL WITH 

MATTER 

To understand the underlying dynamics of the model with matter (i.e., with Q ^ 0) we 
shall shall study the model with two minimally coupled scalar fields together with matter 
having energy density p with the barotropic equation of state given by ( |2.1|) . This model is 
obtained by setting n = 2 in (|2.4|) , whence we obtain the five- dimensional dynamical system 
given by: 

f^ = fi(2g-37 + 2) 

CLT 

— ^ = ^1 g - 2 - ^A;i$i2 
dr 2 
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dr 2 

where 



= *.(«- 2) -^fa*.^ 

dr 2 

— ^ = + l + — /Ci^i) 

ar Z 

"^^^ ^$2(5 + 1 + ^^2*2) (4.1) 



g = + 2*? + 2*^ - - 



and 



= -1 + O + ^1 + ^'2 + $1 + $ 
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A. Invariant Sets and Monotonic Functions 

The zero-curvature models constitute a four-dimensional invariant set. The models with 
no matter also constitute a four-dimensional invariant set. 
The function 

W = - (A 2) 

has derivative 

dW 

^ = 2(2 - 37) W (4.3) 

We observe that this function is monotonic when Q 7^ (i.e., non-zero matter) and (Jl -|- 
^l + '^l + ^l + ^l — l) 7^ (i.e., non-zero curvature). We also observe that the sign of 37 — 2 
signifigantly changes the dynamics of these models. For example, in the case of interest here 
37 — 2 > 0, whence is a decreasing function of time r. This immediately implies that 

• There exist no periodic or recurrent orbits in the full five-dimensional phase space 
(this does not preclude the existence of closed orbits in the invariant sets Q = and 
= 0; however, we shall be primarily concerned with the dynamics of the models in 
the complete phase-space with matter and non-zero curvature). 
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• The future asymptotic state lies within the invariant set f2 = 0. Matter becomes 
dynamically unimportant to the future. 

• The past asymptotic state lies within the set of zero-curvature models. 

B. Stability of Equilibria 

The equilibrium points can be classisfied into two sets; those with J7 = and those with 
Q 0. All equilibrium points listed in Table H exist in this case with = 0, and Table 
m lists the equilibrium points with = together with the additional eigenvalue due to 
the addition of matter. Using the function W above, we can further conclude that those 



equilibrium points in the set i7 7^ necessarily must have zero-curvature. Table |T| lists the 
eigenvalues found in the invariant set 7^ 0. 

Let us focus on the stability of the attractors in the full physical phase space. All late- 
time attractors (sinks) occur in the invariant set = 0. In the previous section we found 
that A and CS are the only sinks in the invariant set = (clearly, all of the saddles 
remain saddles in the full five-dimensional phase space). The additional eigenvalue for the 
equilibrium point A in the full physical phase space is given in Table |I| and is negative if 
Yl'i=i k^'^ > 3^. But this is always satisfied when Yl'i=i ^-i"^ > \ and 7 > |, and hence A is 
a sink and assisted inflation is a global attractor. Similarly, from Table |I| the equilibrium 
point CS IS always a sink for Y^i=i < \ and hence the two-field curvature scaling solution 
remains the global attractor in this case. 

The early-time attractors lie in the zero-curvature invariant set and consist of massless 
scalar field models. From Table |I| we see that the massless scalar field models corresponding 
to the repelling equilibrium points MSF are always sources (for 7 < 2). 
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C. Matter Scaling Solutions 



In the case of a single scalar field there exist zero- curvature FRW 'matter scaling' so- 
lutions when the exponential potential is too steep to drive infiation, in which the scalar 
field energy density tracks that of the perfect fiuid (so that at late times neither field is 
neglible) In it was shown that whenever these matter scaling solutions exist they 
are the unique late-time attractors within the class of fiat FRW models . The cosmological 
consequences of these scaling models have been further studied in For example, in 

these models the scalar field energy density tracks that of the perfect fiuid and a significant 
fraction of the current energy density of the Universe may be contained in the homogeneous 
scalar field whose dynamical effects mimic cold dark matter; the tightest constraint on these 
cosmological models comes from primordial nucleosynthesis bounds on any such relic density 
PJ9|J2^. The stability of these flat, isotropic matter scaling solutions was studied within the 
class of spatially homogeneous cosmological models with a barotropic perfect fluid and a 



scalar fleld with an exponential potential in |[T0|. It was found that while the matter scaling 
solutions are stable to shear perturbations, for realistic matter with 7 > 1 they are unstable 
to curvature perturbations. 

Returning to the models under investigation here, none of the equilibrium points with 
Q ^ can be late-time attractors for 7 > |. Indeed, from Table |T| all such equilibrium 
points are seen to be saddles. In particular, the two-fleld matter scaling solution correspond- 
ing to the equilibrium point MS, which exists for XlLi ^ ^ saddle. From Table 



nil we see that the flrst eigenvalue associated with MS is positive, while the real parts of 



the remaining four eigenvlaues are all negative. This is consistent with the stability analysis 
of matter scaling solutions in models with a single scalar fleld which found that the models 



were unstable to curvature perturbations when 7 > | [|T0[- However, these two-fleld mat- 
ter scaling solutions may still be of physical import. We note that when the curvature is 
zero, the two-fleld matter scaling solution is an attractor (all four eigenvalues of MS in the 
four-dimensional zero-curvature invariant set have negative real parts - so that MS is a sink 
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in this invariant set), as in the case for the matter scahng solution in a single field model. 
Note also from Table |ITT| that both of the single-field matter scaling solutions, corresponding 
to the equilibrium points MSi and MS2, have two positive eigenvalues, so that again the 
solution with multiple scalar fields is the 'stronger' attractor. 



V. QUALITATIVE ANALYSIS OF THREE-SCALAR FIELD MODEL 

Let us now consider models with more than two scalar fields. For simplicity, we shall 
exclude a matter term here. However, from the previous section we easily determine the 
essential properties resulting from the inclusion of a matter field. In particular, in this case 
a monotonic function exists and this enables us to prove the qualitative results outlined 
below. Let us begin with the three-scalar-field model, obtained by setting n = 3 and i7 = 
in (|2.4|) . In this case the resulting six- dimensional dynamical system is given by: 

T / ON ^7 ^ 2 

— — = -q/Jq - 2) ^2*2 

— ^ = ^3 g - 2 - 4-A;2$3' 
dr 2 

— ^ = $1 g + l + — A;i^i 
dr 2 

—± = $2(g+l + _A;2^2) 

-i = $3(g+l + ^A;2Vl/3) (5.1) 
dr 2 



where 



and 



q = 2^l + 2^1 + 2^^ - - $^ - $^ 



6H^ 

Again it would possible to choose simplified variables as in pB[ via a rotation in field 



space as was done in recent work |IJJ21||. However, we shall not do this here. Indeed, we 
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shall not present a complete qualitative analysis similar to that done in section |T|, since the 
essential features are similar and the detailed analysis would be long and painful. Rather, 
let us describe the main effects of including a third scalar field on the inflationary solutions. 

There exists a zero-curvature assisted inflationary solution which now corresponds to the 
equilibrium point given by 

where 

In this solution all of the three scalar fields scale together at late times. The corresponding 
eigenvalues are 



-2, 



2 ' 
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- (^{K^ - 6) ± ^{K^-Qy + 8K^{K^-6) 
^ ({K^ - 6) ± ^{K^ - 6)2 + %K\K^ - 6)) . 



It is known |jT3[ to be a stable late-time attractor for all parameter values for which the 

solution is inflationary (i.e., < 2; recall the point does not exist if > 6). 

There are three solutions in which two scalar fields scale together asymptotically and 

the third is negligible. Assuming that the third scalar field is zero (\l/3 = $3 = 0), the 

coordinates of the corresponding equilibrium point, denoted by Puo, are given by 

r hh' k,^k2 ^6(fc,2 + k,') - k,%' ^Qjk,' + k,') - k,%' 

^ VQik,^ + k,'y VQik,^ + k^^y ' ' VQik.' + k^') ' ' VQik,' + k^') ' ^' 

(5.2) 

Four of the eigenvalues are given by ( |3.5| ), which all have negative real parts. 

There are three solutions in which one scalar field scale dominates asymptotically and the 
remaining two are negligible. Assuming that the first scalar field is non-zero (\E'i ^ ^ $1), 
the coordinates of the corresponding equilibrium point, denoted by Pioo, are given by 



Two of the eigenvalues are negative, one is positive, and there is an additional eigenvalue 
which is negative ii kf < 2 and positive ii 2 < kf < 6. 

In both of these cases the additional (remaining) two eigenvlaues can be calculated and 
are given by 

{g-2<0,g + l>0}, 

where q is the deceleration parameter evaluated at the equilibrium point. Hence, the point 
Pi2o is a saddle with one eigenvalue with positive real part. The equilibrium points denoted 
by Pio3 and P023 are also saddles with one eigenvalue with positive real part. In additional, 
the point Pioo is a saddle with two eigenvalues with positive real parts (if kf <2, and three 
eigenvalues with positive real parts \i kf > 2). The same is true for the equilibrium points 
denoted by P020 and Poo3- 

Consequently there is a nested set of equilibrium points. At the top is the stable three- 
scalar field assisted inflationary solution. In the next layer there are three two-scalar fleld 
models which are saddles with one eigenvalue with positive real part. In the final layer there 
there are three one-scalar field models which are saddles with two eigenvalues with positive 
real parts (or three eigenvalues with positive real parts). Associated with this dynamical 
nesting are cosmological models with very interesting physical properties. 

This will follow through in the case of n scalar fields. There will be a unique stable 
n-scalar field assisted infiationary solution. There will then be n of the (n — l)-scalar field 
models which are saddles with one eigenvalue with positive real part. There will be |n(n — 1) 
of the {n — 2)-scalar field models which are saddles with two eigenvalues with positive real 
parts. And so on. Finally, there will be n of the (l)-scalar field models which are saddles 
with n — 1 (or n — 2) eigenvalues with positive real parts. As one 'goes up' the nested 
structure the equilibrium points respectively become 'stronger attractors' (i.e., the stable 
manifold of the equilibrium points increases in dimension) . 

There is also a three-field curvature scaling solution corresponding to the equilibrium 
point given by 
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V]/ ^ $ ^ 

whose associated eigenvalues are given by 

-1 ± v/l + 4K-2(2-/s:2), -1 ± VSz, -1 ± VSi 

This equihbrium point is a sink whenever > 2, in which case it represents an FRW model 
with negative curvature (else it is a saddle and represents a positive curvature model) . 

Finally, there are saddle equilibrium points corresponding to the Milne model and 
the one- and two-field curvature scaling solutions, and a set of equilibrium points with 
{Xli=i = 1) = 0} corresponding to massless scalar field models, a subset of which are 
sources. 

A complete qualitative analysis can be done for n scalar field models. All of these results 
can be proven by induction (see, for example, [|14|)- The n-scalar field assisted infiationary 
solution is given by |]13 



R(t) oc tP 

and 

ki4>i = kj(f)j;\fl < i j < n, 

and 

n ^ 

We note that in the two-scalar field model, although infiation can occur for potentials that 
are steeper than in the single-field case, it cannot occur for arbitrarily steep potentials. For 
example, if ki = k2 = k, then infiation occurs if k'^ < 4. However, for n-fields, if ki = k for 
all i, then infiation occurs if fc^ < 2n] e.g., k'^ < 8 for four scalar field models. 
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VI. CONCLUSION 



We have studied multi-scalar field FRW cosmological models with exponential potentials, 
extending previous analysis by including non-zero curvature and barotropic matter. We 
have used dynamical systems techniques, and by establishing a monotonic function in the 
complete dynamical phase space (which includes both matter and curvature), we have been 
able to deduce global results. 

In section |T| a comprehensive qualitative analysis was presented in the case of two scalar 
fields with no matter. We concluded that the two-field assisted inflationary solution A is 
the global attractor when Yl'i=i ^i"^ > \ the two-field curvature scaling solution CS 
is the global attractor when X]i=i < |- ^ subset of the massless scalar field solutions 
MSF are always the early-time attractors. Consequently, we found that in all cases both 
scalar fields are non- negligible in generic late-time behaviour; this is an interesting and 
unexpected result and is contrary to the commonly held belief that in multi-field models 
with exponential potentials the scalar field with the shallowest potential would dominate 
at late times (indeed, we have shown that the single field power-law inflationary models 
always correspond to saddles). We note that both the assisted inflationary solution and the 
massless scalar field early-time attractors correspond to zero-curvature models. However, the 
curvature is not always dynamically negligible asyptotically because the two-field curvature 
scaling solution has non-zero curvature. 

The zero-curvature assisted inflationary FRW scaling solutions [jTSf are of particular 
importance since, through the combined effect of multiple uncoupled scalar fields, each 
having an exponential potential, power-law inflation is possible even when each individual 
scalar field need not be a source for inflation. We have discussed the stability of the two-field 
assisted inflationary model, and generalized previous results by including non-zero curvature 
to show that for an appropriate range of parameter values the assisted inflationary solution 
is the global late-time attractor. For these parameter values the single field power-law 
solutions Pi and P2 were shown to be saddles, and we showed that there are allowable 
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parameter values for which either Pi and P2 are both inflationary, or one is inflationary 
while the other is not, or and both are non- inflationary, perhaps leading to new interesting 
physical scenarios. 

In section ^ we studied the two-scalar field model with barotropic matter. A monotonic 
function was established in the resulting phase space. This proved that the matter must 
be negligible at late times and we found that A and CS are the only global sinks and 
that consequently assisted inflation and the two-field curvature scaling solution are the 
global late-time attractors in their appropriate respective parameter ranges. This confirmed 
the result that both scalar fields must be dynamically non-negligible in generic late-time 
behaviour, and establishes the stability of the two-field assisted inflationary model when 
matter is included. The monotonic function also shows that the early-time attractors lie in 
the zero-curvature invariant set, and we showed that they consist of a subset of the massless 
scalar field models. 

For 7 > |, all of the equilibrium points with Q ^ were shown to be saddles (see Table 
pn] ). The two-field matter scaling solution corresponding to the equilibrium point MS was 
shown to have a single positive eigenvalue. Both of the single-field matter scaling solutions, 
corresponding to the equilibrium points MSi and MS2, were shown to have two positive 
eigenvalues, so that again the solution with multiple scalar fields is the 'stronger' attractor. 
We note that when the curvature is zero, the two-field matter scaling solution is the late- 



time attractor, consistent with the stability analysis in []I0[. These matter scaling solutions, 
and particularly the two-field matter scaling solutions, give rise to new transient dynamical 
behaviour and may be of physical import. For example, there are solutions which spend 
a period of time with the scalar field mimicking the barotropic fluid in which there is a 
non-negligible scalar field (dark matter) energy density (corresponding to a matter scaling 
saddle equilibrium point) and subsequently evolve towards a scalar-field dominated power- 
law inflationary epoch (corresponding to a single-field saddle equilibrium point or a two-field 
assisted inflationary attractor) with an accelerated expansion, perhaps explaining current 
high redshift data. 
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In section we discussed three- and multi-scalar field models (where, for simplicity, a 
matter term was excluded). In the three-scalar field model we again established the assisted 
inflationary solution and three-field curvature scaling solutions as the stable late-time at- 
tractors. We then considered n-scalar field models, and established a nested structure for 
the m-field scaling (assisted inflationary) solutions. The n-scalar field assisted inflationary 
solution is again the late-time attractor. All of the m-field (with m < n) scaling solutions are 
saddles and in general the equilibrium points corresponding to the m-field scaling solutions 
will have n — m eigenvalues with positive real parts so that the equilibrium points corre- 
sponding to the greater number of non-negligible scalar fields are, respectively, the 'stronger 
attractors'. Again we should emphasize that Malik and Wands showed that the multi- 
scalar field assisted inflationary solution is a late time attractor by utilizing a rotation in 
field space; indeed, the stable modes in a general stability analysis of this solution are pre- 
sumably related to the isocurvature perturbations orthogonal to the attractor trajectory in 
field space obtained in their analysis. 

Finally, from previous investigations |Q] of spatially homogeneous scalar field cosmological 
models with an exponential potential and barotropic matter and from the above analysis, we 
can conclude that the assisted inflationary solution is a global attractor for all ever-expanding 
spatially homogeneous multi-field cosmological models with exponential potentials provided 
Yll=i^i'^ > 1/2- We can also conclude that the multi-field curvature scaling solution is 



a 



global attractor for models of Bianchi types V and Vllh provided ^^^=1^1 < 1/2 p9 



and a multi-field generalization of the Feinstein-Ibanez anisotropic single-field solution ||30 
is the global attractor for models of Bianchi types III and Vlh if X^ILi ^ V^- Indeed, 
there will be n-field generalizations corresponding to all equilibrium points of the single-field 
Bianchi Type B models (cf. [Q). 

Finally, we note that spatially flat FRW matter scaling solutions also exist in the context 
of generalized assisted inflation. In [^] it was shown that in the higher- dimensional context, 
in the six-dimensional model the assisted dynamics between the scalar fields mimics the 
behaviour of a relativistic fluid (7 = 4/3), while for higher dimensions the scalar fields 
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dominate the radiation component, perhaps leading to a 'moduh' problem for the early 
universe. 
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APPENDIX: DETERMINATION OF THE EIGENVALUES FOR THE 

EQUILIBRIUM POINT A 



If we let X = (\E'i, \E'2, $i, $2) then we are able to write the dynamical system ( |3.1|) as 
(iX/(ir = F(X) where F is an analytic function from R'^ R^. Standard results from 
dynamical systems theory states that the local behaviour near an equilibrium point, Xq, 
of a system of non-linear autonomous differential equations of the form c/X/rfr = F(X) is 
equivalent to the linearized system dX^/dr = DF(Xo)X in a neighborhood of the equilibrium 
point provided the eigenvalues of the matrix DF(Xo) have non-zero real parts [^. For the 
system given by equation (|3.1|) the derivative matrix has the form 



DF(X] 



g-2 + 4^i2 
4\Efi$2 



4-^1 \E'2 



4^2 $1 



$2 {4:^2 + — k2) 



-2 *2$1 

g-2$i2 + l + ^A;i^i 



-2^1 $2 
-2 ^2'^'2 - VGk2 $2 



-2$i<l>2 



-2 $1 $2 

As an example, we shall evaluate DF(X.) at the assisted inflationary equilibrium point ( p^ 
which is denoted as A. 

2 



g-2<l>2' + l + ^fc2^^ 



DF{A) 



2 K' 



3 k,^ 
2 



-3 



3 ki /c2 
1 K{-AK^ + 3ki^] 
6 ^7 

2 K^K 



3 ki k2 

1 2 K'^ 

2 3 k2^ 
2 K^K 



1 K{K^-3k,^] 
3 



1 K^K 



ki' 
1 K^K 



3 ki ^2 
where K = y/R"^ 



3 ki k2 

1 K {-A + 3 k2'^) 
6 h? 



3 ki k2 

"3 

1 ^2(6-^2) 
3 ki k2 



3 ki k2 

1 K{K^-3k2''] 
3 fc? 

3 ki k2 

"3 k? 
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The characteristic polynomial of the matrix DF{A) can be simplified to the form 

The roots of ca{X) yield the eigenvalues seen in equation (3.5). The eigenvalues of all the 
other equilibrium points can be computed in a similar fashion as above. 
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TABLES 

TABLE L Equilibrium points in the two-scalar field model with no matter. In the table, A and 
CS correspond to the two-field assisted inflationary solution and the two-field curvature scaling, 
solution respectively. 



Solution/Label 


Coordinates 


Deceleration 


Curvature 


Eigenvalues 




{<I'l,*2,*l,*2} 


Parameter, q 







Assisted Inflation, A 1'2, *2} (see qa (see 3.3) 



(see 3.5) 



Power-Law, Pi 



{-^,0,,/l-^,0} 



fc/ - 2, 



2 2 



Power-Law, P2 



k2 I A,- 2 
{0,-^,0, 



^2^-2 



2 r, '^2 '^2 



2 2 



, 2 2 , 2(fc 2 -I- fc 2) _ ^ 2t 2 

Curvature Scaling, CS" {- — ,- — , , ^ ^ ,9,9 — liVSi, 



3ki 3k2 VSki V3k2 



k-^k^ 



-1±-^1 + A[2{k-^ + k-^)-l\) 



Curvature Scaling, CSi {- — ,0, ^^,0} 

3fci vSfcl 



2- fci2 



-2, 1,-1 ± Jl + Ak~^{2-k^) 



Curvature Scaling, CS2 {0, - — , 0, — ^} 

3fc2 V3fc2 

Milne, M {0,0,0,0} 



2 k2 
^2 



-2, 1, -1 ± W 1 + 4fe2'^(2 - k^) 



-2, -2, 1, 1 



Massless Scalar Field, MSF {-^10,6 J 1 - 1'q2,o, 0} 



where < 'I'q2 < i 



0, 4, 3 + — fci-^o, 3 H ek2\J 1 - 



where e = ±1. 
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TABLE II. Equilibrium points with Q = in the two-scalar field model with matter. Each 
equilibrium point has = and the coordinates given in Table |. The additional fifth eigenvalue 
is displayed. 

Label A Pi P2 CS CSi CS2 M MSF 

5*'' Eigenvalue T^^t^ " 37 ^'i^ ^ 37 fcj^ - 87 2- 37 2- 37 2- 37 2- 37 3(2 - 7) 



TABLE III. Equilibrium points with 7^ in the two-scalar field model with matter. Note 
that in each case '^R = and q = (87 — 2)/2. 



Solution/Label 


Coordinates 


Eigenvalues 




{r2,1'i,*2,'I'i,*2} 





FRW, F {1,0,0,0,0} 37-2, |(7-2),f{7-2),|7,f7 

37 v/67 „ 767(2-7) „i ^,^_2)L 

((7 - 2) ± ^(7 - 2)2 + 87(7 - 2)[1 - 37fcr']) 

37 , V67 , ^67(2 - 7) 1 ^'7-2) 

(^(7 - 2) ± ^(7 - 2)2 + 87(7 - 2)[1 - 37^2^']) 



Matter Scaling, A/Si {1 ir,-- — ^,0, -^^ — ^,0} -(7 - 2), -7, 37 - 2 

^ ^ k^^ 2ki 2ki ' 2^' ^' 2 ' ' 

3 ' 
4 



Matter Scaling, MS2 {1 " ' 0' " ' 0' "Mt; 77(7 - 2), -7, 37 - 2 

fej 2fc2 2A;2 2 2 

3 



4 



Matter Scaling, MS {1 - 37 {k{''^ + fcj"^) ^ _ ^ 37 - 2, - ("{7 - 2) ± \/(7 - 2)2 + 87(7 - 2) 

V / 2A:i 2k2 4 \ • 

\/67(2 - 7) 767(2 - 7) 1 3/, ^, _^ /; , ^ , — , -2 



2fci 2fci 



} J (^(7 - 2) ± 7(7 - 2)2 + 87(7 - 2)[1 - 37(fci-' + k 



2 
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